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I. INTRODUCTION 

An accurate knowledge of the top-quark mass is im- 
portant. Combined with other standard model param- 
eters, the top-quark mass can be used to constrain the 
expected standard model Higgs mass. The most precise 
constraints place an upper bound on the Higgs mass of 
154 GeV/c 2 [l| at the 95% confidence level. Further, 
since the top quark will be produced in copious quanti- 
ties at the LHC, its mass can serve as a benchmark for 
6-jet energy calibrations, which are otherwise difficult to 
study in data. 

Historically, top-quark mass measurements have been 
limited by uncertainties in the simulation of jet energy 
measurements. For example, in the CDF Run I measure- 
ment [2] the top-quark mass systematic uncertainty was 
4.9 GeV/c 2 , of which the contribution from the uncer- 
tainty on jet energy measurements was 4.4 GeV/c 2 . In 
comparison, for a modern CDF Run II measurement [3| , 
the systematic uncertainty is 1.2 GeV/c 2 , of which the 
contribution from jet energy uncertainties is 0.7 GeV/c 2 . 
This dramatic improvement in the jet energy uncertainty 
is made possible by exploiting the hadronic IF-boson de- 
cay. By constraining the reconstructed IF-boson mass 
to agree between data and simulation, the average sim- 
ulated jet energy bias is determined and applied to all 
jets in the event. In this technique, widely used in Run 
II top-quark mass measurements, most of the jet energy 
uncertainty becomes a statistical uncertainty rather than 
a systematic uncertainty. The remaining systematic un- 
certainty results from the assumption that simulation is 
biased by the same amount, on average, for the 6-jets as it 
is for the VF-boson jets that are used for the calibration. 

While jet energy uncertainties are no longer as large as 
they once were, they still represent a significant fraction 
of the total uncertainty on the top-quark mass. Further, 
this dramatic improvement in the world average mass res- 
olution depends upon the reliability of a single technique. 
Thus, it is desirable to develop independent methods to 
measure the top-quark mass as a cross-check. A novel 
technique has been proposed for measuring the top-quark 
mass in a manner that is almost completely independent 
of the calorimeter using the measured transverse distance 
4] that fe-hadrons from the top quarks travel before de- 
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caying (Lxy) The transverse momenta of the top- 
quark decay products depend approximately linearly on 
the top-quark mass. In turn, this means that the lifetime 
and decay length of the 6-hadrons depend approximately 
linearly on the top-quark mass. When this measurement 
was performed using 695 pb _1 of integrated luminosity 
[f|, it became the first top-quark mass measurement to 
be mostly independent of calorimeter-based uncertain- 
ties, however it was also very statistically limited. Be- 
sides tripling the integrated luminosity with respect to 
the previous measurement, we also improve the statisti- 
cal resolution by incorporating more information about 
the event. Similarly to 6-hadrons, the transverse mo- 
menta (pt) of leptons from VF-boson decays also depend 
linearly on the top-quark mass, and are mostly indepen- 
dent of the calorimeter jet energy scale. Since the mo- 
mentum of the leptons is mostly uncorrelated to that of 
the 6-quarks, this is complementary information, and it 
is an ideal variable to add to the measurement, as has 
been previously proposed Q 

In this paper we present measurements of the top- 
quark mass using both the lepton transverse momentum 
and the mean decay length variables. Both measure- 
ments are performed upon the same events which pass an 
event selection that is designed to isolate ti events where 
the W-boson from one top-quark decays to two quark 
jets, and the other decays leptonically to an electron or 
a muon plus a neutrino (the lepton plus jets channel). 
This decay channel was chosen because it occurs more 
often than any other final state that can be isolated with 
a high signal purity. We compute the mean Lxy of the 
leading two 6-tagged jets and the mean transverse mo- 
mentum of the leptons identified in the events. We mea- 
sure the top-quark mass through comparisons with the 
mean Lxy and mean lepton px from analysis simulations 
( "pseudoexperiments" ) performed for a variety of top- 
quark mass hypotheses. To illustrate the dependence of 
these variables on the top-quark mass, the expected dis- 
tributions for ti signal events passing our event selection 
are shown in Figure [TJ It should be noted that our mea- 
surements are sensitive to very different event character- 
istics than typical mass measurements, and thus require 
unique treatments. In particular, it is more important for 
us to correctly model the boost of the top quarks than 
for other measurements. Thus, we reweight our simu- 
lated signal events to a more accurate parton distribu- 
tion function. Further, the Lxy values that are measured 
in our simulation are sensitive to a variety of possible 
modeling inaccuracies. We determine correction factors 
to be applied to the simulated Lxy values by comparing 
the Lxy measurements between data and simulation for 
bb events, parameterized as a function of the jet energies. 

Our measurements will be mostly independent of 
calorimeter-based jet energy uncertainties. However, 
some small dependence will remain because a loose jet 
energy cut is applied when selecting events, as will be 
explained below. Further, this parameterization of the 
Lxy correction as a function of jet energy introduces a jet 
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FIG. 1: Distributions of our measurement variables for simulated tt events after passing our event selection for top-quark 
masses separated by 50 GeV/c 2 . These distributions are normalized to unit area. 



energy uncertainty to the measurement. In order to mini- 
mize this latter effect, we develop and apply an algorithm 
to measure jet energies for our Lxy correction parame- 
terization that is based upon charged particle tracking. 

This paper is organized as follows. We first present 
the relevant parts of the CDF detector in Section Hi] and 
the event selection used in this analysis in Section IIIII 
We then explain the procedures for determining the nor- 
malizations as well as the shapes of the Lxy and lepton 
Pt distributions for our backgrounds. In Section [V] we 
explain how we calibrate our signal distributions. In Sec- 
tion I VII we present the procedures we apply to extract 
the top-quark mass from our final Lxy and lepton mo- 
mentum distributions and present our results. Finally, 
in Section IVIII we explain in detail how our systematic 
uncertainties are estimated. We conclude with some pro- 
jections of the future potential for this type of measure- 
ment. 



II. THE CDF II DETECTOR 

We describe the most important parts of the CDF de- 
tector for this analysis. A more complete description 
can be found elsewhere Q. The CDF detector consists 
of cylindrically symmetric layers of hardware, each de- 
signed to perform specific functions. The CDF tracking 
system consists of a silicon tracker, inside of a wire drift 
chamber, immersed in a 1.4 T magnetic field produced 
by a superconducting solenoid. Calorimeters encircle the 
tracking system and consist of alternating layers of scin- 
tillator and absorber. The calorimeter system is split into 
an inner calorimeter that is designed to absorb and mea- 
sure the energies of photons and electrons, and an outer 
calorimeter that is designed to absorb and measure the 
energies of hadrons. 

The wire drift chamber consists of 96 wire planes 
grouped into eight superlayers covering the radial region 
between 0.40 and 1.37 m from the beam axis, and less 



than 1.0 in pseudorapidity. The superlayers alternate 
between axial and ±2° stereo angle to provide both axial 
and longitudinal hit information. The inner tracker con- 
sists of eight layers of double sided silicon, covering the 
radial range between 1.35 and 29 cm, and less than 2.0 
in pseudorapidity. Again, to provide longitudinal hit in- 
formation, some layers contain double-sided silicon that 
combines axially oriented strips with strips oriented ei- 
ther at 90° or 1.2° relative to the axial direction. The 
silicon tracker is vital for vertexing, providing a two- 
dimensional impact parameter resolution of about 40 /jm 
for isolated tracks with high transverse momentum, in- 
cluding uncertainties on both the track trajectory and the 
location of the collision vertex. The wire tracker is vital 
for the lepton measurements in this analysis, providing 
a transverse momentum resolution of about 5% for 50 
GeV/c muons. Leptons and &-jets are reliably identified 
out to a pseudorapidity of about 1.0, after which the ef- 
ficiency drops rapidly due to tracks falling outside of the 
range of full tracker coverage. Muon candidates are iden- 
tified by matching tracks in the inner tracking detector 
with track stubs in the outer muon drift chambers. Elec- 
tron candidates are identified by matching tracks in the 
inner tracking detector with showers in the electromag- 
netic calorimeter. An energy resolution of about 3% is 
achieved for electrons with 50 GeV of transverse energy. 



The Tevatron proton-antiproton bunches cross one an- 
other at a rate of 2.5 MHz. During each crossing one 
or more collisions are likely to occur. In order to reduce 
this flow of information to a manageable level and to 
select collision events of interest to particular analyses, 
CDF employs a three level triggering system to sequen- 
tially reject events that are less relevant. The final event 
stream is read out at a rate of roughly 100 Hz. Details 
of the triggering system are provided in Q. 
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III. EVENT SELECTION 

In comparison to the previous publication using only 
Lxy [|| , we tighten the event selection in the analyses pre- 
sented here to reduce the systematic uncertainties. The 
statistical sensitivities of both the lepton momentum and 
the decay length measurements depend linearly on the 
fraction of events in which top quarks are produced 
and scale as the square root of the number of events se- 
lected. This tightened selection improves the former and 
worsens the latter, and has little impact on the final sta- 
tistical sensitivity. 

The data used in this analysis were collected between 
March 2002 and May 2007, and correspond to an inte- 
grated luminosity of 1.9 fb _1 . Events passing the full 
trigger and event selections described below are stud- 
ied to determine the expected event counts and un- 
certainties for each signal and background type. Top 
quarks almost always decay to Wb. Our selection crite- 
ria are designed to accept events where one VF-boson 
decays to an electron or muon plus neutrino and the 
other decays to two jets. We start from a triggering 
stream that requires one electron (muon) to have trans- 
verse energy (momentum) greater than 18 GeV (GeV/c); 
these requirements will later be tightened slightly. Once 
events are accepted by the trigger, they are saved, re- 
constructed, and studied in greater detail. Calorime- 
ter towers are clustered together within a cone of radius 
R = ^((f>tow - 4*jet) 2 + (vtow - Vjet) 2 = 0.4 [l(| to form 
jets. At least three jets must be found with |?7| < 2.0, 
and transverse energy greater than 20 GeV after correct- 
ing for multiple interactions, calorimeter response, noise, 
and other non-uniformities. No attempt is made to cor- 
rect for interactions from spectator partons ( "underlying 
event") or out-of-cone effects [l(J, but systematic uncer- 
tainties are assigned for these effects. To account for the 
neutrino and suppress the QCD background, a quantity 
called the missing transverse energy is used, which is de- 
fined as the transverse component of the four-momentum 
vector that is needed to conserve momentum in the event 
[Tl| . The missing transverse energy in the event must be 
greater than 20 GeV. Additionally, an electron (muon) 
must be identified with transverse energy (momentum) 
greater than 20 GeV (GeV/c). Electrons are identified 
from a track pointing at a cluster in the calorimeter which 
matches the expected shape profile. Most of the energy 
of this cluster is required to be confined in the electro- 
magnetic calorimeter, the track momentum is required 
to agree with the measured calorimeter energy to within 
a factor of two, and if the track is consistent with the 
electron having originated from a photon conversion, the 
electron candidate is vetoed. Muons are formed from 
tracks in the muon chambers which are matched to tracks 
found in the inner tracking system. The calorimeter en- 
ergy deposits along the muon trajectory must be consis- 
tent with that of a minimum ionizing particle. Calorime- 
ter isolation is based upon the fraction of the lepton's 
energy (fi SO ) in a cone of radius R=0.4 centered on the 



lepton, excluding the energy in the calorimeters from 
the lepton itself. Both electrons and muons must sat- 
isfy fiso < 0.1. This cut eliminates most jets that are 
mistakenly identified as leptons ("fake leptons") as well 
as real leptons which result from 6-hadron decays. Fur- 
ther, at least one collision vertex must be reconstructed 
from tracks in the event, and the track of the lepton must 
pass within 3 cm along the beam axis of the highest mo- 
mentum vertex to minimize contamination from multiple 
interactions and tracking errors. One electron or muon is 
required to pass these cuts, but to suppress events from 
the dilepton channel the event is vetoed if any other lep- 
tons are found passing a much looser set of cuts which 
also allow non-isolated and forward leptons. Finally, one 
or more jets must be identified (tagged) as originating 
from a 6-quark, as explained below. In the case where 
only three jets pass our selection, two of them must be 
tagged in order to reduce the large VF+jets and non-IF 
QCD backgrounds in this kinematic region. 

Jets containing 6-quarks are identified using a tagging 
algorithm known as SecVtx (l2| . which relies upon the 
long lifetime of hadrons originating from 6-quarks. This 
algorithm attempts to construct a secondary vertex us- 
ing tracks that are likely to have originated from a b- 
hadron decay. If a vertex can be found which is signifi- 
cantly displaced from the primary vertex, then the jet is 
tagged. Specifically, the SecVtx algorithm selects tracks 
associated with the jet that are displaced from the fit- 
ted primary collision vertex, and that are well resolved 
in both the silicon and the outer tracking chamber. As 
a first pass the algorithm uses relatively less stringent 
cuts on track selection, and attempts to fit at least three 
tracks into a displaced vertex. If no vertex is found the 
tracking requirements are tightened significantly and the 
algorithm searches for a two-track vertex. Further cuts 
are applied to eliminate tracks that are likely to originate 
from material interactions or long lived strange particles. 
If a secondary vertex is found the jet is tagged as a 6-jet 
if the transverse distance from the primary vertex pro- 
jected onto the jet direction (Lxy) and its uncertainty (a) 
satisfy Lxy/c > 7.5. It should be noted that charmed 
daughters of the 6-hadron are also likely to travel a signif- 
icant distance before decaying. The SecVtx algorithm is 
deliberately designed to be loose enough to attach some 
tracks from these tertiary decays into one "pseudovertex" 
at a position that is averaged between two real vertices. 
Since the boosts of the charm hadrons depend on the 
boost of the 6-hadron, this extra information does not 
dilute our mass resolution. 



IV. SAMPLE COMPOSITION 

For this analysis we will need to know the Lxy and 
lepton pt distributions for each of our signal and back- 
ground samples, as well as their relative normalizations 
after full event selection is enforced. Specifically, we 
normalize our signal and background distributions us- 
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ing the results of a cross section measurement that was 
performed on the same data using identical event selec- 
tion. The procedures for this cross section measurement 
are briefly described in subsection IIVAI The procedures 
that are used to determine our Lxy and lepton px shapes 
for the backgrounds are deferred to Section HVBI 



A. Sample normalization 

In this section we give a brief outline of how the cross 
section measurement is performed to determine our sig- 
nal and background normalizations along with the asso- 
ciated uncertainties. Further information can be found 
in the publications of similar cross section measurements 
[HI [3 , or in the Ph.D. thesis [l5[ which explains the 
technique we apply in detail. For this measurement tech- 
nique we begin by determining the expected backgrounds 
assuming the expected standard model ti cross section. 
The ti cross section will then be revised and the back- 
grounds redetermined in an iterative manner until the 
number of expected events matches the observed event 
counts both before and after 6-tagging. 

First we determine the numbers of events for some of 
the rarest processes from simulation. Backgrounds mod- 
eled from simulation include single top production as well 
as the electroweak diboson (WW, WZ, ZZ) and Z plus 
jets final states. The diboson events are simulated using 
PYTHIA version 6.216 [l6[. The single top and Z plus 
jets samples are simulated using other programs (MadE- 
vent [13] for single top, and ALPGEN version 2.10 prime 
[lH for Z plus jets). Hadronic showers are simulated 
using pythia. For each sample, the event count is de- 
termined and scaled according to the theoretical cross 
section, branching ratio, and detector and trigger accep- 
tances. The ti signal is also simulated using pythia. Its 
cross section is initially set equal to its standard model 
expectation. For all of these samples, EvtGen is used 
to determine the lifetimes and masses of the various 6- 
hadron species within jets. A further scaling is applied 
to correct the 6-tagging efficiency modeling based upon 
data-driven studies |14| . 

Since fake leptons are difficult to simulate accurately, 
the non-IL QCD contribution is determined from data. 
Missing transverse energy templates are made for ti, W 
plus jet, and QCD distributions, and used in a fit to de- 
termine the QCD normalization. Templates for the fake 
electrons are filled from events where isolated electron 
candidates are selected and required to have a calorimeter 
shower profile that is consistent with QCD fakes rather 
than true electrons. For muons the standard cuts are 
kept except that the isolation cut is inverted to require 
greater than 20% isolation instead of less than 10%. Dif- 
ferent binnings, fit ranges, and cuts are applied and the 
differences in the results are taken as a systematic un- 
certainty. This fit is done separately in the tagged and 
pretagged cases. 

The remainder of the observed pretagged events are 



all taken to originate from W plus jets. The W plus 
jets sample is simulated using ALPGEN version 2.10 
prime [18|, and the resulting partons are showered us- 
ing pythia. This sample is simulated in various bins of 
jet multiplicity for the Wbb, Wcc, Wc, and W plus light 
flavor final states. These samples are then weighted by 
their theoretical cross sections and combined. It is impor- 
tant to make sure that the proper number of heavy fla- 
vor jets are simulated before 6-tagging is applied. Pythia 
showering will double count heavy flavor production from 
the ALPGEN simulation, so this overlap must be removed. 
Further, since the samples are only generated at leading 
order, a correction is needed to rcweight the heavy flavor 
fraction. This correction is taken from comparisons be- 
tween data and simulation in related samples with higher 
statistics 

By construction, this procedure yields pretagged back- 
ground and signal normalizations that exactly match the 
observed number of data events. After 6-tagging, how- 
ever, fewer events were predicted than were observed, 
and this discrepancy is attributed to events. Thus, the 
signal cross section was increased and the analysis was 
repeated in an iterative fashion. The ti cross section at 
which we found agreement was 8. 2 ±0.7 pb (ignoring lumi- 
nosity uncertainties). The resulting event counts passing 
full event selection and tagging requirements are shown 
in Table HI along with the number of jets in these events 
which were 6-tagged and included in our Lxy analysis. 



TABLE I: Estimated signal and background contributions af- 
ter full event selection. 



Source 


Events 


Recorded b- Tags 


Wbb 


25.4 ± 7.0 


37.5 ±9.7 


Wcc or Wc 


13.9 ±4.6 


15.9 ±4.7 


VLplus light flavor 


16.9 ±3.7 


17.9 ±3.7 


non-WQCD 


18.8 ± 12.7 


20.2 ± 13.2 


Electroweak 


9.0 ±0.4 


11.3 ±0.5 


Single Top 


8.4 ±0.4 


13.4 ±0.7 


ti 


478.3 ± 40.3 


659.3 ±45.5 


Total 


570.8 ± 44.3 


775.5 ± 50.2 



B. Background Lxy and lepton pr shapes 

The dominant backgrounds for our analysis are W plus 
heavy flavor (6 and c jets), W plus light flavor mistags, 
and non-VL QCD events. Along with single top and di- 
boson events, these distributions and the signal account 
for about 99% of events passing selection. The remaining 
events come from the Z plus jets background, for which 
the related W plus jets background Lxy and lepton pt 
distributions were used. 

Single top samples (with masses rat = 165, 170, 175, 
and 180 GeV/c 2 ) were generated in MadEvent and de- 
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cayed in pythia, in both the s- and t- channels. Results 
for the s- and t- channels were combined, weighted ac- 
cording to their expected theoretical cross sections, and 
the final decay length and lepton momentum distribu- 
tions were fit to Gaussian plus exponential distributions. 
The trends in the fit parameters were extrapolated to 
other mass points, and were used to generate new Lxy 
and lepton px distributions for each mass hypothesis. 
The W plus jets background is selected in exactly the 
same way as for the cross section measurement, as is the 
background for non-W QCD electrons. For non-VF QCD 
background in the muon channel, instead of using non- 
isolated muons, we selected muon candidates with high 
energy deposition in the calorimeter, or that were asso- 
ciated with tracks that were highly displaced from the 
collision vertex. These cuts were chosen based on studies 
suggesting that they produced a minimum of bias in the 
lepton pt distribution. 

Our final Lxy and lepton px distributions are shown 
in Figure [2] for events passing full selection. As cross- 
checks the same distributions in the background domi- 
nated one and two jet bins are shown in Figures [3] and 
HI These cross-check samples are used in evaluating the 
background based systematic uncertainty as described in 
Section El 



addition, since the CTEQ5L PDF overestimates the frac- 
tion of the proton momentum carried by colliding quarks, 
the ti events produced within the quark-antiquark anni- 
hilation channel have artificially high boosts. To com- 
pensate for this effect, different parton distribution func- 
tions were studied. The Les Houches Accord PDF In- 
terface [22j was used to determine the probabilities for 
each parton involved in the collision to have the gener- 
ated momentum for a given PDF. One can then reweight 
the simulated events to construct distributions appropri- 
ate to a different PDF as explained in (23|. We use this 
prescription to reweight all of our ti events to the next-to- 
leading order CTEQ6M [24| parton distribution function 
for gluon fusion and quark annihilation events separately. 
A further weighting is then applied to gluon fusion events 
to scale the gluon fusion fraction to the value expected for 
the sample's top-quark mass (20% for m t = 150 GeV/c 2 , 
10% for m t = 200 GeV/c 2 ). These combined reweight- 
ings result in new distributions to be used in our mass 
measurement and lead to a 1.7 GeV/c 2 shift in the top- 
quark mass for the lepton pt analysis and a 0.9 GeV/c 2 
shift for the Lxy analysis, relative to the values obtained 
using CTEQ5L. A similar prescription is used to reweight 
to other PDFs and gluon fractions to evaluate systematic 
uncertainties, as will be explained in Section IVlTl 



V. CORRECTIONS TO THE SIGNAL 
SIMULATION 

Our ti events are simulated in pythia for various hy- 
pothetical top-quark mass values. A number of correc- 
tions are applied to the simulated results as discussed 
below. In Section IV Al we will discuss corrections that 
are needed to account for simulated parton distribution 
function inaccuracies. It is also necessary to correct the 
simulated Lxy measurements to match observations from 
data. An overview for this procedure is given in Sec- 
tion [VB] We will parameterize the Lxy correction as a 
function of the jet transverse energy as measured using 
tracking. The procedure for this will be discussed in Sec- 
tion IV CI We will summarize the complete decay length 
calibration procedure in Section IV Dl 



A. Parton distribution functions 

For this analysis it is vital to have accurate models of 
the lepton and jet boosts. Since these quantities depend 
on the energy of the colliding partons, an accurate model- 
ing of parton distribution functions (PDFs) for particles 
within the proton is very important. The ti samples were 
generated using the leading order CTEQ5L parton distri- 
bution function 20]. One drawback of this PDF is that it 
underestimates the rate of ti production by gluon fusion 
(5% in simulation observed versus about 15% expected 
[211]). Since gluon fusion produces events with slightly 
smaller boosts on average, this will bias the analysis. In 



B. Lxy calibration strategy 

A number of effects may bias the decay length mea- 
surement in simulation. Inaccuracies in the EvtGen [l9| 
values for the hadron lifetime or in the simulated pro- 
duction fractions would have a direct impact. Similarly, 
inaccuracies in the PYTHIA fragmentation model would 
lead to the wrong boost, and thus the wrong average de- 
cay length, of the 6-hadrons. Entirely different problems 
may arise from any inaccuracies in the modeling of the 
tracking system, which could lead to biases in the vertex- 
ing results. Our approach is to calibrate the simulation 
directly to the data to compensate for all of these biases 
simultaneously. Systematic uncertainties on this calibra- 
tion will be discussed in Section IVII Fl We select our 
calibration sample so that the tagged jets will be almost 
exclusively 6-jets. We select dijet events where the jets 
are required to open at a wide angle with A(j) > 2.0. Both 
jets must be 6-tagged, and one of the jets must contain 
a well resolved muon with at least 9.5 GeV/c transverse 
momentum. We also apply additional cuts to minimize 
overlap between jets which will be discussed later. We 
estimate these samples to be about 95% bb, with the re- 
maining 5% coming from charm contamination, which is 
accounted for as a small systematic uncertainty. 

Since jets from bb events tend to have a much lower 
transverse energy than those from ti events, the possibil- 
ity that a different calibration is needed for higher energy 
jets must be taken into account. To this end, we bin our 
bb jets according to their energy, and derive the needed 
correction bin-by-bin, which we apply based on the mea- 
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FIG. 2: Signal, background, and data for the Lxy and lepton pt distributions passing full event selection under hypothesized 
top-quark masses that are close to the measured results. The left plot is for the Lxy measurement, using top-quark mass 
168 GeV/c 2 , and the right plot is for the lepton pr measurement, using top-quark mass 173 GeV/c 2 . 




FIG. 3: Background prediction compared with data (black points) in the one-jet control region for Lxy (left) and lepton pr 
(right). 



sured energy of the tagged signal jet. Great care must 
be taken in doing this, however, as this kinematic based 
correction directly introduces a jet energy scale uncer- 
tainty to the analysis. If, for example, the simulation 
underestimates the energies of jets, then the average de- 
cay length in a given energy bin will be too high in the 
simulation, throwing off the calibration. This is an un- 
avoidable uncertainty for any kind of calibration using 
dijets. In fact, even if we could convince ourselves that 
it was unnecessary to parameterize this decay length cal- 
ibration as a function of energy, a significant jet energy 
uncertainty would still be needed to cover the determi- 
nation of which jets pass selection in the simulation. As 
our goal is to minimize the calorimeter jet energy uncer- 



tainties, we choose to parameterize our calibration based 
upon the energies of jets measured using tracking rather 
than the calorimetry. 



C. Track-based jet energies 

We develop a straightforward algorithm for computing 
the track-based energy of a jet. We start with one of 
our jets that has been clustered in the calorimeter in the 
usual manner. We select tracks that are within R = 0.4 
of the calorimeter jet direction in r\ — (f> space. The tracks 
themselves are required to be well resolved in both the 
wire tracking chamber and in the silicon, and they must 
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have a transverse momentum of at least 1 GeV/c. We 
also require that the tracks pass within 3 cm along the 
beam axis of the fitted primary collision vertex. This cut 
eliminates about 85% of the contamination from multiple 
interactions. We then take the total transverse compo- 
nent of the sum of our track four- vectors as our tracking 
transverse energy. We make no attempt to correct for the 
missing neutral particles. Our goal is not to measure the 
true jet energy, but rather to measure a quantity that is 
proportional to the true energy (in this case, the charged 
transverse energy) in a manner that agrees well between 
data and simulation. 

We calibrate the energies of our track jets with a simi- 
lar approach to the one used for calorimeter jets at CDF 
(icj |. We select events where one photon and one jet 
are found back-to-back (A<fi > 3 radians), where no ex- 
tra jets in the event above 3 GeV in Et are allowed, 
and strict cuts on photon quality are applied to mini- 
mize fakes. Under such a selection, the transverse mo- 
mentum of the photon should reflect the true transverse 
momentum of the jet. As a first step we consider the ra- 
tio of the track-based transverse energy of the jet to the 
measured transverse energy of the photon. The distri- 
bution of this ratio shows good agreement between data 
and simulation as seen in Figure 5(a) for a particular 
range of photon energies, and this agreement also holds 
for higher energy photons. The mean fraction of the pho- 
ton transverse energy that is carried in the track-based jet 
transverse energy is shown for a range of true jet trans- 
verse energies (determined by the energy of the photon) 
in Figure |5(b)| The extent of the agreement between 
data and simulation is given by the ratio of these trends 
and is shown in black in the same figure. A line is fitted 
to these ratios and represents the calibration that will 
be applied to the track jet energies that we measure in 
our simulation. Specifically, for a given tagged jet in our 



simulated bb or ti sample, we begin by measuring the 
corrected transverse energy of the jet in the calorimeter 
and its associated track-based jet transverse energy. We 
then correct this track-based transverse energy accord- 
ing to its calorimeter-based transverse energy and the 
fitted function in Figure [5 (b)| Uncertainties in the simu- 
lation of the calorimeter-based jet energy measurements 
are found to have a negligibly small impact on the correc- 
tion factor that is determined in this manner. There are, 
however, significant statistical uncertainties on the fitted 
function that translate into a systematic uncertainty on 
our results as will be explained in Section IViTl 

When applying this calibration to the simulated bb and 
ti samples, however, it is important to account for the 
fact that these events contain more jets than those in the 
photon calibration sample, and tracks from other jets 
may fall into the jet cones and bias the measured track 
jet energies upwards. If our decay length calibration pa- 
rameterization proves to have a trend over jet energy, 
then such biases in the ti sample will have a direct im- 
pact on the correction factor that is applied to a given 
jet. On the other hand, in the bb samples there will only 
be a bias to the extent that pythia does not properly 
model the amount of jet overlap that is observed in data. 
We start by minimizing this problem as much as possible 
by removing events where a second jet is close to the jet 
we are studying and might contribute overlapping parti- 
cles in both data and simulation. For our bb events we 
veto events where either of our tagged jets is between 
AR of 0.7 and 1.2 of any other jet with energy greater 
than 9 GeV. This cut was chosen to remove most events 
that might have extra jets in a region that could overlap 
our primary jet without eliminating jets with energetic 
out-of-cone QCD radiation. 

In the ti simulation, we do not veto events in this 
manner. Instead, we develop a correction procedure to 
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FIG. 5: (a): The distributions of track-based jet transverse energy divided by the photon pr for events with photon pr between 
30 GeV/c and 40 GeV/c. (b): The mean fraction of track-based jet transverse energy / photon pr for Pythia (red) and the 
data (blue) as a function of the measured photon transverse momentum (assumed to be the true jet transverse energy) . The 
ratio of these trends between data and simulation is fitted to a line which is then used to correct measured track jet energies 
in the simulation. 



remove the effects of tracks from other jets that over- 
lap our jet cones. We do this by plotting the track mo- 
mentum associated with 6-tagged tt jets as a function of 
the AR between the track and the jet direction as mea- 
sured in the calorimeter. We plot distributions binned in 
the generator-level px for the parent 6-hadron in our jet. 
The higher the energy of the 6-hadron, the narrower the 
AR distribution, as expected. We fit these distributions 
to two components, a "primary" part, and an "overlap" 
part. The tracks originating from the 6-hadron and asso- 
ciated fragmentation products (the primary jet) are mod- 
eled by a Gaussian multiplied by a Fermi function to force 
the momentum fraction to converge to zero at small AR. 
The contribution from underlying event, minimum bias, 
and other jets (overlap), turns out to be well modeled by 
a quadratic function in AR, multiplied by a Fermi func- 
tion to account for damping effects as tracks pass out of 
the active detector range. The fits are performed sepa- 
rately for each region of 6-hadron px- Examples of these 
fits (along with cross-check fits for bb events) are shown 
in Figure [HI We then use these fit results to extract cor- 
rection factors, parameterized by the 6-hadron's px, to 
remove the average amount of momentum from charged 
particles expected to fall inside the track jet cone from 
other sources. 

There are a number of approximations and assump- 
tions that have gone into these vetoes and corrections. 
For example, the fitting functions could be inappropri- 
ate, or the vetoes applied to the bb events might be inad- 
equate. If so, then using an alternate cone size to select 
track jets will lead to a direct systematic shift. Thus, the 
systematics for this procedure are evaluated by changing 
the track jet cone size to 0.7, removing overlap accord- 
ing to the revised jet size, and repeating the analysis, as 



explained in Section IVII Fl However, two other cross- 
checks are also performed to improve confidence in the 
procedure. First, we have no physical motivation for us- 
ing a quadratic function as the base shape of our over- 
lap, so in one cross-check we repeat the analysis, mod- 
eling our overlap as though it were distributed perfectly 
uniformly in r\ — 4> space (a line times a Fermi function 
in AR). There is, of course, no physical basis for us- 
ing this symmetric distribution either, since complicated 
correlations between jets according to tt kinematics and 
sculpting from the jet clustering algorithm could render 
the shape asymmetric. This alternate shape is simply 
a cross-check that leads to an overestimated amount of 
overlap falling inside of the jet cone. When the analysis is 
repeated with the alternate shape it leads to a top-quark 
mass measurement result that is shifted by an amount 
that is smaller than the systematic we will eventually 
assign for the overlap removal procedure. As a second 
cross-check we look at the tagged jets that are selected 
to be back-to-back with the muon jet in the bb sample. 
If the overlap fitting function is incorrect, then we would 
expect the shape of the primary part of the fit to be bet- 
ter modeled in these bb events, since the overlap is about 
a factor of four smaller than in tt events. These bb com- 
parison distributions are shown as the dashed curves in 
Figure [5] By far the worst agreement is seen in the low- 
est bin of hadron px (shown on the left). In every other 
bin of 6-hadron momentum, the differences between the 
bb and tt results can barely be distinguished, as is the 
case in the right hand plot. But since less than 8% of tt 
jets have hadron momenta that fall into the lowest bin, 
the systematic mass shifts caused by enforcing the alter- 
nate bb shapes in the primary distribution will be well 
less than the systematic that we will end up claiming. 
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FIG. 6: These plots show the charged jet momentum spread in AR for pythia ti events. Specifically, the AR of each track 
relative to the calorimeter jet direction is plotted as the black points for tagged ti jets, weighted by the track momentum. The 
results are then fit to determine the fractions of primary jet energy and overlap energy from other jets for ti. The solid blue 
curves are the fit shapes for the primary and overlap fit components and the solid red curve is the full fit result. The dashed 
blue curve shows the expected primary distribution based upon fits in the bb sample, (a): results for b-hadron Pt less than 
20 GeV/c. This is the only kinematic range where there is any significant disagreement between the ti and bb results. Much 
more typical is (b): results for fe-hadron Pt between 50 GeV/c and 60 GeV/c. 



D. Final decay length calibration 

With the track jets in hand we can evaluate the de- 
cay length calibration parameterization in the bb sample 
as described above. The calorimeter is still used to se- 
lect jets and determine their direction, but we loosen the 
calorimeter jet energy cut significantly, and exclude most 
jets with energies near this cut by applying a track jet 
transverse energy cut. As a result, calorimeter driven 
jet energy uncertainties will play a minimal role. Ad- 
ditionally, we apply the overlap vetoes and corrections 
as described above. Since it is possible for a jet to pass 
our tracking energy cuts while failing the calorimeter en- 
ergy selection, fluctuations within standard, calorimeter- 
based jet energy uncertainties can still cause events to 
pass in and out of selection. However this is a small ef- 
fect which is only significant in the lowest energy track 
jet bins, and thus leads to a small systematic uncertainty. 
Figure [7] shows the trends in the mean decay length for 
data and simulation. The ratio of these trends gives us 
the correction that should be applied to the measured 
decay lengths in data, depending on the measured track 
jet energy in the signal samples. The distribution of the 
track jet energies to which the calibration will be applied 
is overlaid. 
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FIG. 7: The decay length calibration parameterization versus 
track jet transverse energy is shown here. It is developed from 
the solid black data and simulated average decay lengths, the 
ratio of which gives the Scale Factor points, which are then fit 
to determine the calibration. Also shown is the distribution 
of jets to which the calibration is applied. For the Lxy data 
and simulation the vertical scale represents the mean decay 
length in cm, while for their Scale Factor ratio and the tt 
distribution this axis is unitless. 



VI. MASS MEASUREMENT METHOD 

Given our mean measured Lxy and lepton pt values 
from the data (associated with the distributions of Fig- 



ure [2]) we need to determine the associated top-quark 
masses and statistical uncertainties. We simulate exper- 
iments under a variety of hypothetical top-quark masses 
and use them to perform each of the measurements as 



13 



described below. 



B. Measurement using both variables 



A. Single variable measurements 

Pseudoexperiment events are drawn from the events 
that are used to construct the signal and background 
distributions where the probability of each event is given 
by its PDF weighting as discussed in Section [Vj The 
mean Lxy (px) of the tagged jets (leptons) in these events 
will be used to measure the mass. Samples of tt events are 
generated under 23 hypothetical top-quark mass values 
ranging from 140 GeV/c 2 to 220 GeV/c 2 , and the decay 
length results are corrected according to the track jet 
energies by the fit results of Figure [7] 

Uncertainties from the background normalization are 
small and are wrapped into the pseudoexperiments. A 
total of 92.5 ±17.1 background events are expected, how- 
ever due to the iterative nature of the tt cross section 
evaluation, there is a small additional uncertainty on the 
background normalization of ±3.7 events based upon the 
theoretical uncertainty on the input tt cross section, lead- 
ing to a total background count uncertainty of ±17.5. For 
each pseudoexperiment, the total number of background 
events is fluctuated according to a Gaussian with the 
above mean and RMS to determine an expected back- 
ground normalization. The resulting number is then 
fluctuated according to Poisson statistics to determine 
a background normalization for each pseudoexperiment. 
Given a fixed number of observed data events, the excess 
is taken to be signal. Further, for each pseudoexperiment 
the Lxy calibration parameterization is fluctuated within 
the fitted statistical uncertainties shown in Figure [7J 

The mean Lxy and lepton px pseudoexperiment re- 
sults are Gaussian in shape and are fit to a Gaussian for 
each hypothetical top-quark mass. Examples of these fits 
are shown in Figure [5) To evaluate the top-quark mass 
results for the Lxy and lepton pt measurements, the cen- 
tral values of these Gaussians are plotted as a function 
of top-quark mass and are fit to a quadratic polynomial. 
The mean Lxy and lepton px values measured in data are 
then converted to the measured top-quark mass values 
according to this polynomial. To extract the statistical 
uncertainties, the central Lxy and lepton px pseudoex- 
periment values are shifted up and down by the standard 
deviation of these Gaussian fits. Then the difference be- 
tween the measured top-quark mass according to the un- 
shifted polynomial and the top-quark masses resulting 
from these one standard deviation shifts are taken to be 
the asymmetric one sigma statistical uncertainties on the 
measurements. The fitted polynomials are shown in Fig- 
ure OH The mean Lxy and lepton px values measured in 
data are also shown as the horizontal black lines, along 
with the projections that are used to determine the sta- 
tistical uncertainties. 



The pseudoexperiments from the single variable results 
are used to plot two-dimensional mean Lxy versus mean 
lepton px distributions. The results for the two most 



extreme mass hypotheses are overlaid in Figure 10(a) 

The observed data produce a point on this two- 
dimensional plane. Given this point, our task is to de- 
termine the most likely value of the top-quark mass and 
the associated statistical errors. To accomplish this, we 
evaluate a likelihood for each mass hypothesis according 
to the data. The likelihood is simply the probability that 
if the true mass were the one in our hypothesis, the mean 
Lxy and lepton px results would fluctuate as far away or 
farther than the results we see in the data point. This 
probability is taken from pseudoexperiment results such 
as those shown in Figure 10(a) Specifically, we evaluate 
a "distance" that our data point is from the expected 
central value for our mass hypothesis, and take the like- 
lihood that the hypothetical mass is correct to be the 
fraction of pseudoexperiments which are "farther away" 
from the expected values than our data point. 

For this approach to be meaningful, a reasonable def- 
inition of "distance" must be used. We choose our def- 
inition so that distances are equal for points along the 
equal probability contours of a two-dimensional Gaus- 
sian centered at the expected mean Lxy and lepton px 
values, with the expected standard deviations. Here, the 
expected means and standard deviations are taken from 
the Gaussian fits to the mean Lxy and lepton px pseu- 
doexperiment results that were described for the single 
variable measurements. Then, the "distance" from the 
expected value is defined according to the metric in Equa- 
tion [TJ 
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Here, 8P t (5L xy ) is the difference between the mean 
lepton px (Lxy) of the data and the fitted central value 
for the hypothesized top-quark mass, and <jp t {&L X y) is 
the fitted standard deviation of the mean lepton px (Lxy) 
for the hypothesized top-quark mass. In principle this 
approach could be modified to account for a correlation 
between the two variables, but this is unnecessary as the 
correlations are empirically determined to be negligibly 
small. 

Based on this equation, the fraction of pseudoexperi- 
ments for each hypothesis for which the distance metric 
is evaluated to be larger than that for the data is taken 
as the likelihood for the hypothesized mass. Finally, the 
likelihood values for each mass point are plotted with sta- 
tistical uncertainties determined by adding in quadrature 
the statistical uncertainties from the pseudoexperiments, 
and the statistical uncertainty due to the number of sim- 
ulated events from which they are drawn. We run enough 
pseudoexperiments (4000) that the size of our simulated 
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FIG. 8: Pseudoexperiment distributions and fit results for example hypothetical top-quark mass values, (a): Mean Lxy 
pseudoexperiments for mt = 165 GeV/c 2 . (b): Mean lepton pr pseudoexperiments for mt = 173 GeV/c 2 . The results of these 
fits will be used to construct Figure [9] 
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FIG. 9: The mean values of the Gaussian fits to the pseudoexperiment results from Figure|S]are fit to the quadratic polynomials 
plotted here in red. The one sigma statistical uncertainties from these fits lead to the blue contours. The mean values from data 
are shown as the horizontal black lines. These values are then translated into top-quark masses according to their intersections 
with the red polynomial, and into asymmetric statistical uncertainties according to their intersections with the blue polynomials. 



tt samples is the primary limitation. These likelihoods 
are fit to a Gaussian as shown in Figure [T0(b)| The mean 
of the Gaussian is taken to be the result of our combined 
measurement with a statistical uncertainty given by the 
RMS of the fit. 

C. Results 

Data events passing event selection have a mean Lxy 
of 0.590 ± 0.017 cm and a mean lepton pr of 55.2 ± 
1.3 GeV/c. Based upon these values, the mass measure- 
ment with the decay length technique yields a result of 



166.9^8 5 GeV/c 2 , and with the lepton transverse mo- 
mentum technique yields a result of 173. 51^9 GeV/c 2 , 
where the errors are statistical only. For the simultane- 
ous measurement with both variables the fit to data is 
shown in Figure [10] and corresponds to a mass result of 
170.7 ±6.3 GeV /c 2 . 

Some sanity checks were run for each of the three top- 
quark mass measurements. For these checks, nineteen 
additional top-quark mass samples were generated with 
top-quark masses varying between 152 and 193 GeV/c 2 . 
Ten of these samples were blind (the masses were hidden 
from the authors until after the measurements were fin- 
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FIG. 10: (a): distribution of mean Lxy versus mean lepton pr from pseudoexperiments for extreme mass cases of 140 GeV 
and 220 GeV. (b): the pseudoexperiment results for the 23 mass points considered are used to determine the likelihood of 
agreement with the data according to the metric of Equation [1] and the results are plotted and fitted here. The mean fit result 
is taken as the measurement result, and the RMS represents our statistical uncertainty. 



ished). Pseudoexperiments were thrown using these sam- 
ples, and the means of the measured top-quark mass re- 
sults proved to be consistent with expectations to within 
the statistical uncertainties, indicating the method is un- 
biased. Further, the pulls (defined as the difference be- 
tween each measured mass and the generated mass di- 
vided by the statistical uncertainty) of the pseudoexperi- 
ment results prove to have a width that is consistent with 
1.0, indicating that the estimated statistical uncertainties 
are reliable. 



VII. SYSTEMATIC UNCERTAINTIES 

Next we will discuss the systematic uncertainties for 
this measurement, and the motivations for the proce- 
dures we apply. A list of our final systematic uncertain- 
ties can be found in Table [VTl 



A. Background uncertainty 

As mentioned above, the uncertainty on the back- 
ground composition and shape is evaluated in the control 
regions of the one- and two-jet bins, shown in Figures [3] 
and HI The differences between the observed and ex- 
pected means are shown in Table [TTJ The largest dis- 
agreements are observed in the one-jet bin of Lxy, and 
in the two-jet bin for lepton px- Since both of these 
worst-case shifts are larger than their uncertainties, they 
are taken as the uncertainties on the background mean 
results and are scaled by the background fraction to de- 
termine the systematic errors. 



TABLE II: Background shifts and uncertainties in the one 
and two-jet control regions. The uncertainties account for 
both statistical effects due to data limitations as well us un- 
certainties in the relative contributions from the individual 
backgrounds from the cross section measurement. 

Variable Shift 

One- Jet Lxy -0.0131 ± 0.0082(cm) 

One- Jet lepton p T 0.25 ± 1.25(GeV/c) 

Two- Jet Lxy 0.0022 ± 0.0118(cm) 

Two- Jet lepton p T -2.22 ± 1.05(GeV/c) 



B. QCD radiation uncertainty 

Uncertainties in the simulation of QCD radiation could 
have a significant impact on the results of an analysis 
like this which is dependent upon an accurate model- 
ing of the boost of the decay products of the ti system. 
Comparisons of the dilepton boost for Drell-Yan events 
have been made between data and simulation [23[ and 
used to constrain inaccuracies in the modeling of initial 
state radiation in quark-antiquark interactions, pythia 
parameters were varied to conservatively bracket the pos- 
sible disagreement between data and simulated initial 
state radiation, and analogous parameters were simulta- 
neously varied for the final state radiation by an equiv- 
alent amount. Signal events were generated with these 
parameters shifted up and down, and these samples were 
compared with each other and the nominal pythia sam- 
ple. Half of the largest mass shifts between any pair of 
these three samples was taken as the QCD radiation sys- 
tematic uncertainty. 
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C. Parton Distribution Function uncertainty 

As described in SectionfVl for this analysis we reweight 
our events to match both the predictions of the next to 
leading order CTEQ6M parton distribution function, and 
the expected gluon fusion top production fractions pre- 
dicted by theory. The advantage of using the CTEQ6M 
PDF is that it includes a prescription for estimating PDF 
uncertainties [24J- The degrees of freedom of the PDF 
can be parameterized in 20 orthogonal sources of un- 
certainty that are commonly called eigenvectors. For 
each of these uncertainty sources there are two parton 
distribution functions, where the parameters associated 
with these eigenvectors are shifted up or down to cover 
a 90% confidence interval. We reweight our top-quark 
mass 175 GeV/c 2 sample to each of these forty alter- 
nate PDFs, taking half the full mass shift for each pair 
as a systematic uncertainty, and adding them in quadra- 
ture. While these uncertainties are intended to represent 
a 90% confidence interval, we conservatively use them as 
one sigma systematics instead. While we fix the fraction 
of ti events produced by gluon fusion interactions to theo- 
retical expectations, it should be noted that some of these 
eigenvector variations are expected to change this frac- 
tion. Thus, we allow the gluon fractions to float around 
their expectations for purposes of determining systematic 
uncertainties on the PDF results. 

One uncertainty that is not accounted for by these 
eigenvectors is the uncertainty on the strong coupling 
constant, a s (m z ) = 0.1176 ± 0.0020 [25j]. To study the 
effects of this uncertainty we reweight to the CTEQ6A 
and CTEQ6B PDFs [26] , which are two different series of 
PDFs constructed with varying a s values in intervals of 
0.002. We average the mass shifts obtained when varying 
the PDFs, and arrive at an uncertainty that is roughly 
half as large as the eigenvector uncertainty. We observe 
consistency between the A and B series PDFs. We add 
this uncertainty in quadrature to the eigenvector uncer- 
tainty to determine our full PDF systematic. As a final 
cross-check to the results of the CTEQ Collaboration, we 
also reweight to the MRST Collaboration's NLO PDF 
MRST2004 [27[. We observe agreement well within our 
stated eigenvector uncertainty when this result is com- 
pared with the corresponding CTEQ6A/B PDFs. 

D. Generator uncertainty 

In CDF top-quark mass analyses it is conventional to 
re-evaluate the top-quark mass u sing samples produced 
with the herwig 6.510 generator [2g|, and take the shift 
from the PYTHIA mass result as a systematic uncertainty. 
Note that many of the differences between these gener- 
ators will double count our existing systematic uncer- 
tainties. The different fragmentation models between 
PYTHIA and herwig will double count the decay length 
scale, jet energy, and QCD Radiation uncertainties, her- 
wig also does not properly handle QED radiation off 



of leptons from the VF-boson decays, which is instead 
inserted with the PHOTOS [2^| program. Differences 
in these approaches will double count our lepton energy 
scale uncertainty. The generators also have minor differ- 
ences in the applied top width (and herwig has a sharp 
cutoff preventing the presence of top quarks in the high 
and low mass tails), and only herwig properly handles 
spin correlation between the two top quarks. Despite the 
double counted uncertainties, we follow the convention 
of other analyses for consistency by taking the differ- 
ence between our PYTHIA and herwig mass results as 
a generator systematic. For the Lxy and combined mea- 
surements the statistical uncertainty on our mass shift is 
greater than the shift itself, so we take the uncertainty 
on the shift as our systematic instead. 

E. Lepton momentum uncertainty 

The modeling of the lepton momentum in simulation 
is tested by fitting the invariant mass of Z's in data and 
simulation, separately for electrons and muons. A Brcit- 
Wigner function is used to model the inherent Z width, 
which was convoluted with a Gaussian to account for de- 
tector resolution. Additionally, to model the kinematic 
reduction in the cross section for higher mass Z produc- 
tion this function is multiplied by a decaying exponential. 
When a function modeling a QCD background shape is 
included the fits return zero for its normalization as ex- 
pected due to the high purity achieved by the lepton 
selection. The fit distributions are shown in Figures [Til 
and[U 

The centers of the Breit-Wigner fit results, shown in 
Table IIII1 were compared between simulation and data. 
To evaluate the systematic, the mean lepton pt of the 
signal was scaled by the ratio of the data and simulated 
means for electrons and muons separately, and the shift 
in the measured mass results was taken as a systematic 
(the statistical uncertainties on these fit results are negli- 
gible). Clearly, the disagreements in the electron results 
dominate this uncertainty. 



TABLE III: Centers of the Breit-Wigner functions from the 
fits to the Z peaks shown in Figures [TTJ and [12] 



Sample 


Mass Peak (GeV/c 2 ) 


Muon Simulation 


91.16 


Muon Data 


90.96 


Electron Simulation 


91.81 


Electron Data 


90.84 



F. Decay length related uncertainties 

The procedure for calibrating the decay length mea- 
surements in our signal sample is described in sec- 
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tion IV Bl There are many uncertainties which must be 
considered for this calibration, some due to the modeling 
of 6-jets, and others due to the track jet energy measure- 
ments that are used to parameterize the calibration. 

The decay length calibration has a statistical limita- 
tion due to the data and PYTHIA bb sample sizes. This 
uncertainty is folded directly into the pseudoexperiments 
as explained in section IVI1 but its contribution is quite 
small. There are a number of uncertainties on the photon 
plus jet energy calibration technique. The energy scale 
calibration curve has an associated statistical uncertainty 
which propagates through to a mass uncertainty. A sys- 
tematic uncertainty of 1% is taken on the measured en- 
ergy of the photon in the simulation, which corresponds 
to a 1% uncertainty on the measured track jet trans- 
verse momentum. Finally, as described in [lOj, about 
30% of the photon plus jets sample in data is composed 



of QCD dijet production where one of the QCD jets fakes 
a very clean photon signature through a pion or lambda 
decay. This contamination has been determined to have 
a momentum balance discrepancy compared to the pho- 
ton plus jets signal at the 1% level, and so an additional 
0.3% uncertainty is taken on the track jet momentum. 
Finally, after applying the procedures described in sec- 
tion [V] to minimize (in the case of bb) or correct (in the 
case of tt) for jet overlap and underlying event effects, 
the mass is reevaluated using a cone size of AR = 0.7 
instead of 0.4 for the track jets. The resulting shift in 
the mass results is taken as an additional systematic to 
represent out-of-cone and jet overlap uncertainties. 

Another track jet energy uncertainty arises in connec- 
tion with the simulation of the 6-jets. If the EvtGen 
decay tables do not produce the correct distributions of 
charged particles then this will artificially bias any mea- 
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surements of the tracking energy of the 6-jets. The DEL- 
PHI Collaboration has measured the charged decay mul- 
tiplicity of 6-hadrons [3(3], excluding the decay products 
of long lived light-flavor particles and of excited 6-hadron 
to ground state &-hadron transitions, and determined an 
average of 4.97 ± 0.07. We evaluate this number at gen- 
erator level in our samples using the same exclusions 
and arrive at a mean result of 5.05. This discrepancy 
is very slightly larger than the reported error at DEL- 
PHI, and it cannot be explained by uncertainties in the 
production fractions of different 6-hadron types or on the 
semileptonic decay rate. Under the assumption that ex- 
cess tracks will be distributed randomly in the 6-hadron 
rest frame, this discrepancy should directly translate into 
an equivalent discrepancy on the measured energy of the 
component of the track jet originating from the 6-hadron 
decay. This leads to an additional 1.1% uncertainty on 
the measured track jet energies. 

In addition to jet energy effects, other uncertainties 
are considered in relation to the modeling of the physics 
of the bb sample. It is important to minimize and un- 
derstand any charm contamination in this sample. Our 
studies demonstrate that the muon jets in the sample 
are about 95% likely to be &-jets, and 5% likely to be 
charm jets. But they also suggest that the simulation 
slightly underestimates the number of charm jets. This 
requires a minor correction of the decay length mea- 
sured in this calibration sample, and the charm fraction 
is then fluctuated within its fitted one sigma uncertain- 
ties. The small resulting mass shifts are then taken as a 
systematic uncertainty. Another small uncertainty arises 
from the event selection on the muons in our leptonic 
bb sample. If the simulation does not properly model 
the measurement of the muon momentum then higher 
or lower energy muons (corresponding to higher or lower 
decay length vertices) will pass selection. While our Z 
peak fits above suggest a very accurate modeling of iso- 
lated muons, we conservatively take a 1% uncertainty on 
the muon momentum scale, and evaluate the mass shift 
that results from the new set of events passing selection. 
Finally, it is important to understand the uncertainties 
in pythia's modeling of the 6-quark fragmentation. In 
pythia, the energy carried by the hadron after the frag- 
mentation process is modeled with the Bowler function. 
The DO Collaboration has studied LEP and SLD data 
and determined the pythia tune required to reproduce 
their results [3l[ . Samples of tt events were generated 
according to each of these tunes, and the resulting b- 
hadron energies were found to be about 2% higher than 
under the default pythia tune. As expected, this results 
in a proportionally larger mean decay length of our sig- 
nal &-hadrons. However, since this effect also occurs in 
the bb samples, the effects almost exactly cancel one an- 
other out, illustrating the motivation for our calibration 
procedure. The fragmentation fluctuations do, however, 
produce the following minor fluctuations which are not 
canceled out. When events are reweighted to the alter- 
nate fragmentation distributions, it causes small alter- 



ations to the measured track jet energies and raises the 
muon energy distribution slightly. Accounting for all of 
these effects, the larger of the mass shifts between the 
default pythia sample, and the results after reweighting 
to the SLD or LEP results are taken as a fragmentation 
systematic uncertainty . 

A summary of calibration systematic uncertainties for 
the decay length measurement is shown in Table HVl 



TABLE IV: Calibration based top-quark mass uncertainties 
for the decay length measurement. 



Systematic [GeV/c 2 ] 


Lxy 


Lepton pt 


Simultaneous 


Photon Plus Jet Stats 


0.7 





0.3 


Photon pt 


1.4 





0.6 


Photon Background 


0.4 





0.2 


Track Jet Cone Size 


0.8 





0.3 


Ntrk from 6-hadrons 


1.6 





0.7 


cc Background 


0.2 





0.1 


Semilep Muon Pt 


0.3 





0.1 


Fragmentation 


0.6 





0.3 


Total Calibration 


2.5 





1.1 



G. Multiple interactions uncertainty 

There are two respects in which other interactions dur- 
ing a beam crossing may result in a systematic bias. 
These extra interactions are simulated as overlaid mini- 
mum bias events, however this modeling may not be ac- 
curate, resulting in biased jet energies. These effects are 
described in Section IVlI HI for calorimeter-based jet mea- 
surements, and in Section IVII Fl for the tracking based 
jet measurements. 

However another effect must be taken into account. 
The simulation is tuned to an older dataset correspond- 
ing to an integrated luminosity of 1.2 fb~ . The newer 0.7 
fb _1 of data included in the measurement were collected 
at higher instantaneous luminosities with more interac- 
tions per bunch crossing than the earlier data. To study 
this effect we generated a high luminosity tt sample with 
top-quark mass 175 GeV/c 2 . We then segregated the 
events according to the number of collision vertices that 
are reconstructed in the event (which has been shown 
to be approximately proportional to the luminosity [HI). 
While there is no statistically significant trend over num- 
ber of vertices for the Lxy measurement, there is a sig- 
nificant dependence for the lepton transverse momentum 
measurement as shown in Figure 1131 For electrons, a 
small part of this trend is due to particles from other 
collisions falling into the electron cluster, however the 
primary cause of this effect is the isolation requirement 
for the leptons. Since we require the total calorimeter 
energy found around the lepton to be less than 10% of 
the lepton momentum, low momentum leptons are more 
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likely to fail selection in high luminosity events, as illus- 
trated in Figure IT51 

Like the simulation, we segregate the data based upon 
the number of reconstructed collision vertices. As ex- 
pected, the data have about 15% more reconstructed col- 
lision vertices per event than the standard ti samples. 
The high luminosity simulated events are reweighted to 
reproduce the distribution of the number of reconstructed 
vertices in both the standard ti samples and the data 
in turn. These two reweighted results are equivalent to 
each other within statistics for the Lxy measurement, 
however the lepton transverse momentum is significantly 
higher under the luminosity profile of the data. There 
are insufficient statistics at very high luminosities to reli- 
ably correct for this effect in all of the signal ti samples. 
Instead we take the differences between the associated 
top-quark mass results using the luminosity profiles of 
the data and the simulation for the high luminosity sam- 
ple as a systematic uncertainty for each measurement. 
We emphasize that this uncertainty is due to the sim- 
ple logistics of the luminosity profile that was used in 
the simulation, and is not due to any irreducible physics 
effect. For the decay length measurement the statisti- 
cal uncertainty due to the number of generated ti events 
is larger than the observed systematic shift, and so this 
statistical uncertainty is taken as our systematic instead. 



H. Jet energy uncertainties 

Our jet energy uncertainties can be broken down into 
two categories: those arising from the tracking energy 
measurements which impact the Lxy calibration, and 
those from the calorimeter measurements that are com- 
mon to all our analyses. None of the uncertainties rep- 
resent an uncertainty on the determination of a "true" 
jet energy. Rather, they represent uncertainties in the 
modeling of jet energy measurements in simulation. In 
this section we discuss our evaluation of these uncertain- 
ties and explain why they have minimal correlation with 
the calorimeter-based uncertainties that are claimed by 
other top-quark mass analyses. 

The dominant jet energy uncertainties in this analysis 
arise from the track jets. They are listed in Table HVl as 
Photon Plus Jet Stats, Photon pr, Photon Background, 
Track Jet Cone Size, and N tr k from 6-hadrons. Of these, 
the Photon px (energy bias in the calibration photons) 
and N tr k from 6-hadrons (EvtGen decay multiplicity mis- 
modeling) categories are the largest contributions to the 
jet energy uncertainties. The fourth largest uncertainty 
is due to the limited statistics in the photon plus jets 
data and will not present any difficulty in future high 
statistics analyses. 

The third largest uncertainty is due to the size of the 
cone used to construct our track-based jets. This un- 
certainty may have components from a wide variety of 
physical effects, but minimal correlation to the jet en- 
ergy scale uncertainties of other analyses. The most 



significant component comes from an uncertainty in the 
overlap of particles from other jets falling into the jet 
cone, for which no corresponding uncertainty is claimed 
for calorimeter jets. The only systematic components for 
which there are any correlations to calorimeter-based un- 
certainties are the much smaller underlying event, multi- 
ple interaction uncertainties, and out-of-cone uncertain- 
ties. Most of the multiple interaction contributions are 
vetoed by the z- vertex matching requirement. As for the 
out-of-cone uncertainty, it should be minimal due to the 
application of the photon calibration procedure. It will 
only contribute to the extent in which the simulation 
models out-of-cone effects in 6-jets with a different level 
of accuracy compared to the light flavor jets on which 
the calibration is performed. To summarize, of all our 
track jet energy uncertainties, only a small part of the 
0.8 GeV/c 2 (0.3 GeV/c 2 ) systematic uncertainty on the 
Lxy (combined) measurements that is due to the altered 
cone size could have any correlation to the calorimeter- 
based jet energy scale uncertainties we are claiming, or 
those of other analyses. 

The calorimeter-based uncertainties are split into six 
categories, which are assumed to be independent of one 
another and are described in [l(j. Since these are the 
same categories into which the jet energy scale correc- 
tions are split, these uncertainties are sometimes called 
jet energy scale uncertainties. As for the tracking based 
uncertainties, their impact on the decay length measure- 
ment arises based on which jets pass our event selection 
thresholds. These low energy jets which pass in and out 
of our sample as the jet energies are varied within uncer- 
tainties tend to have a small decay length, and therefore 
bias the average decay length of our sample. Unlike track 
jets, however, for calorimeter jets this effect is present 
in both our bb calibration and our main analysis sam- 
ples, and it largely cancels in the final mass determina- 
tion. Such cancellations are what motivated the choice 
of the Lxy calibration procedure. To evaluate these un- 
certainties, we fluctuate the calorimeter energies of the 
jets within these six categories and reevaluate the miss- 
ing energy of the event, keeping track of which jets and 
events pass selection. We take the resulting mass shifts 
as calorimeter-based systematic uncertainties. 

One concern that arose at this point is that we may 
have over-optimized our procedures to create the fortu- 
itous systematic cancellations described above. This con- 
cern would specifically pertain to our out-of-cone jet en- 
ergy uncertainty. The out-of-cone uncertainties are delib- 
erately fixed to be more conservative than the worst case 
scenario disagreements between PYTHiA-data and her- 
wiG-data out-of-cone comparisons as explained in [Icl |. 
But if the out-of-cone disagreement between simulation 
and data is different for the jets in our bb sample than for 
the jets in our ti sample, then by chance the cancellation 
may lead to an artificially small systematic result. 

To investigate this possibility, we must understand dif- 
ferences between jets in the samples. If jets near the se- 
lection threshold were to have identical properties for the 
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FIG. 13: Effects of luminosity on the mean lepton transverse momentum. These results are evaluated at generator level and 
plotted against the number of reconstructed collision vertices. The higher the luminosity is (as measured by the number of 
vertices), the smaller is the number of low energy leptons that pass the isolation requirement. 



bb and ti samples, then disagreements between data and 
simulation would be identical for the samples and the re- 
sulting systematic cancellation resulting from assuming 
identical out-of-cone uncertainties would be appropriate. 
Fortunately, the differences are small. For our purposes, 
the only relevant differences between the jets in our sam- 
ples are that the bb jets used in our decay length calibra- 
tion are required to contain muons, and that some of the 
tagged jets in our backgrounds are light flavor or charm. 
In all other respects the simulated jets we use near our 
selection threshold are similar. There are two systematic 
cross-checks that we run to address these concerns. The 
results will be shown at the end of this section. 

As explained in [lOj, the out-of-cone jet energy un- 
certainties are parameterized based upon the calorimeter 
energy measurement. Since lower energy jets are broader, 
a lower jet energy corresponds to a larger out-of-cone un- 
certainty. Since the muon's energy is mostly lost for jets 
in the bb sample, it can be argued that we are overes- 
timating the out-of-cone uncertainty for these jets. To 
check the impact this would have, we add the muon's en- 
ergy back in and repeat our systematics evaluation. As 
a second check, we consider the possibility that the mis- 
modeling of out-of-cone effects by the simulation could 
be different for heavy and light flavor jets. We check the 
shifts of the uncertainties which occur when we fluctu- 
ate the size of our out-of-cone uncertainties for charm 
and light flavor jets relative to 6-jets within conservative 
constraints determined by jet shape studies. The fluctu- 
ations from these two cross-checks are taken as residual 
out-of-cone systematic uncertainties. We take our full 
out-of-cone systematic as the quadrature sum of direct 
and residual out-of-cone uncertainties. The results are 
summarized in Table [Vj The systematic uncertainties 
from all effects that have been considered are shown in 



Table EH 



TABLE V: Calorimeter Based Jet Energy Uncertainties. The 
residual uncertainties result from possible inaccuracies in the 
cancellation that occurs for our out-of-cone uncertainty. 



Systematic [GeV/c 2 ] 


Lxy 


Lepton pr Simultaneous 


Eta Dependent 


0.06 


-0.08 


-0.02 


Multiple Interactions 


0.17 


-0.01 


0.07 


Calorimeter Response 


-0.14 


-0.07 


-0.09 


Underlying Event 


0.09 


-0.06 


0.01 


Splash Out 


0.15 


-0.10 


0.02 


Base Out of Cone 


0.18 


-0.28 


-0.06 


Out of Cone Residual Uncertainties 


bb Semileptonic 


0.24 


NA 


0.24 


W plus charm/LF 


0.14 


0.22 


0.30 


Final Out of Cone 


0.33 


0.36 


0.24 


Total Calorimeter JES 


0.44 


0.39 


0.32 



VIII. CONCLUSION 

Using a data sample corresponding to an integrated 
luminosity of 1.9 fb _1 , we measure a top-quark mass of 
m t = 166.9±|;| (stat)±2.9 (syst) GeV/c 2 using the mean 
transverse decay length of 6-jets, m t = 173.5l 89 (stat) ± 
3.8 (syst) GeV/c 2 using the mean transverse momen- 
tum of the leptons from VF-boson decays, and m t — 
170.7±6.3 (stat) ±2.6 (syst) GeV/c 2 using both variables 
simultaneously. To date, these results represent the most 
precise measurement of the top-quark mass using an algo- 
rithm that has minimal dependence on calorimeter-based 
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TABLE VI: Final Systematic Uncertainties. The Lxy Cali- 
bration systematic is the quadrature sum of the systematics 
summarized in Table IIVI The Calorimeter JES systematic 
is the quadrature sum of the systematics summarized in Ta- 
ble El 



Systematic [GeV/c 2 ] 


Lxy 


Lepton pt 


Simultaneous 


Background Shape 


1.0 


2.3 


1.7 


QCD Radiation 


0.5 


1.2 


0.7 


PDF 


0.3 


0.6 


0.5 


Generator 


0.7 


0.9 


0.3 


Lepton pt Scale 





2.3 


1.2 


Lxy Calibration 


2.5 





1.1 


Multiple Interactions 


0.2 


1.2 


0.7 


Calorimeter JES 


0.4 


0.4 


0.3 


Systematics Total 


2.9 


3.8 


2.6 



jet energy uncertainties. Because we have improved the 
systematic uncertainty on this measurement by a factor 
of 3.3 compared to the previous measurement [6| , the pre- 
cision remains limited by statistics and therefore could 
improve significantly by the end of Run II. In contrast, if 
a measurement of this type were performed at the LHC, 
the systematic uncertainties would be the true limitation 
as the statistical uncertainties would be negligible. 

While a prediction of the systematic uncertainties in 
such a measurement at the LHC is beyond the scope of 
this paper, we can predict some of the improvements to 
the systematics that can be made at the Tevatron. If 
the simulated events were regenerated using the proper 



luminosity profile, the multiple interactions uncertainty 
would become negligibly small. Preliminary studies [321 ] 
have also suggested a method for calibrating the lepton 
momentum which should improve our lepton px scale 
systematic by more than a factor of two. If both of these 
analysis improvements were implemented, the systematic 
uncertainties for these measurements at CDF would be 
expected to drop to 3.0 GeV/c 2 for the lepton trans- 
verse momentum measurement, and 2.3 GeV/c 2 using 
both variables simultaneously. It appears that an even 
smaller systematic uncertainty may be attainable by per- 
forming the lepton px analysis in the dilepton channel 
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